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Abstract: The effects of friction on the ion acoustic (IA) wave in fully and partially 
ionized plasmas are studied. In a quasi-neutral electron-ion plasma the friction between 
the two species cancels out exactly and the wave propagates without any damping. If 
the Poisson equation is used instead of the quasi-neutrality, however, the IA wave is 
damped and the damping is dispersive. In a partially ionized plasma, the collisions with 
the neutrals modify the IA wave beyond recognition. For a low density of neutrals the 
mode is damped. Upon increasing the neutral density, the mode becomes first evanescent 
and then reappears for a still larger number of neutrals. A similar behavior is obtained 
by varying the mode wave-length. The explanation for this behavior is given. In an 
inhomogeneous plasma placed in an external magnetic field, and for magnetized electrons 
and un-magnetized ions, the IA mode propagates in any direction and in this case the 
collisions make it growing on the account of the energy stored in the density gradient. 
The growth rate is angle dependent. A comparison with the collision- less kinetic density 
gradient driven IA instability is also given. 



PACS No: 52.35.Fp; 52.30.Ex 



1 Introduction 



Multi-component plasmas comprise different species that, in the presence of waves, may 
be in the state of relative macroscopic motion. In such a situation, friction between 
the species may lead to wave damping (though not always, as we are going to show in 
the forthcoming text). For example, neutrals in a weakly ionized plasma represent a 
barrier for electron and ion motion in a wave field. A similar friction appears in a fully 
ionized plasma when the electron and ion components do not share the same momentum. 
The interaction is described by a friction force Fj = rrijnjVji(vj — v{) in the momentum 
equation for the species j. Momentum conservation implies that for its counterpart I, 
Fi = rrtiriiVijivi — Vj), where mjUjUji = miniUij. If the two species j and I have a large 
mass difference, the friction response of the heavier component is typically omitted as 
negligible in the literature. However, this may yield completely wrong results as we shall 
demonstrate in the forthcoming text using the ion acoustic (IA) mode as an example. 

In the presence of high frequency waves u ^> fij = eB /rrii in a plasma placed in an 
external magnetic field Bq = E>oe z , ions will follow nearly straight lines regardless of the 
direction of the wave-number vector k and the magnetic field vector. For electrons, in view 
of the mass difference, the opposite may hold, uj <C Q e = eB /m e , hence they will behave 
as magnetized and their perpendicular and parallel dynamics will be essentially different 
P]. Ions can behave as un-magnetized in the perturbed state also in case of collisions 
provided that Ui > tti even if at the same time Qi > u, or for short wavelengths A < pi, 
Pi = v T i/Qi, v^i = uTi/rrii. In the case of an inhomogeneous equilibrium, with a density 
gradient perpendicular to the magnetic field vector, in the unperturbed state the ions 
may behave as un-magnetized in case of a low temperature, when their diamagnetic drift 
velocity becomes negligible as compared to electrons [for singly charged ions v*i/v* e = 
Ti/T e , where v*j = nTjn'j Q /(qjBonjo), and n'- = drij/dx denotes the equilibrium density 
gradient]. The same holds in the presence of numerous collisions as above, V{ > Qi, when 
their diamagnetic effects are absent too. 

In all these situations, and neglecting the electron polarization drift (inertia- less limit), 
the wave will still have the basic properties of the IA mode. Within the two-fluid theory 
such a mode in an inhomogeneous plasma [that may be called ion-acoustic-drift (IAD) 
mode] may in fact become growing [T|-[3] in the simultaneous presence of collisions and 
the mentioned equilibrium density gradient perpendicular to B . 
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Within the kinetic theory the mode is also growing in the presence of the same density 
gradient and this even without collisions (due to purely kinetic effects), and the physics 
of the growth rate is similar to the standard drift wave instability [4] . It requires that the 
wave frequency is below the electron diamagnetic drift frequency u* e = v* e k±. 

On the other hand, keeping the electron inertia results in the instability of the lower- 
hybrid-drift (LHD) type [2]"[E]- In some other limits the effects of the same density 
gradient yield growing ion plasma (Langmuir) oscillations jS], or growing electron-acoustic 
oscillations [6]. 

In the present manuscript the friction force effects on the IA wave are discussed, both 
for fully and partially ionized un-magnetized plasmas, and for inhomogeneous plasmas 
with magnetized electrons. The latter implies growing modes within both the fluid and 
kinetic descriptions, and in the manuscript these two instabilities are compared. 

2 IA wave in fully and partially ionized collisional 
plasmas 

The equations used further in this section are the momentum equations for the ions, the 
electrons and the neutral particles, respectively: 

m^i ( 4" + Vi ■ V ) Vi = -eriiVcj) - nTiVrii - m i n i i/ ie ({7 i - v e ) - 171^71^^ - v n ), (1) 



dt 

' d_ 
dt 



m e n e + v e ■ V J v e = en e V(j) - KT e Wn e - m e n e v ei {v e - v t ) - m e n e u en (v e - v n ), (2) 



and 



m n n n +v n ■ V^J v n = -KT n Vn n - m n n n v ni (v n - «<) - m n n n u ne (v n - v e ), (3) 

and the continuity equation 

dn 



+ V ■ (rijVj) = 0, j = e, i, n. (4) 

This set of equations is closed either by using the quasi-neutrality or the Poisson equation. 
The differences between the two cases are discussed below. 

2.1 Friction in electron- ion plasma 

The continuity equation (Tjfl) yields 

vn = una/(kno), v el = un el /(kn ), (5) 
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so that the velocity difference in the friction term v e —Vi = if the quasi-neutrality is used. 
The IA mode propagates without any damping. Hence, the friction force in a fully ionized 
plasma in this limit cancels out exactly even without using the momentum balance. The 
physical reason for this is the assumed exact balance of the perturbed densities: what one 
plasma component loses the other component receives, this is valid at every position in 
the wave and no momentum is lost. 

A typical mistake seen in the literature is to take the friction force term for electrons 
only, in the form m e n e v e iV e . This comes with the excuse of the large mass difference, 
so that the displacement of the much heavier ion fluid, caused by the electron friction is 
neglected. In the case of a fully ionized electron-ion plasma this yields a false damping of 
the IA mode within the quasi-neutrality limit: 



On the other hand, if the Poisson equation is used instead of the quasi-neutrality, one 
obtains [9] 



Here, we have used the momentum conservation v ie = m e v ei lmi and v 1 = c 2 s + 
fde = v Te /u pe . The physical reason for damping in the present case is the fact that the 
detailed balance rin = n e \ does not hold, because of the electric field which takes part 
for small enough wave-lengths. It can easily be seen that for any realistic parameters the 
second term in the real part of the frequency in Eq. (JTj) is much below unity and the 
mode is never evanescent. However, in partially ionized plasmas (see below) this may be 
completely different. 

2.2 Friction and collisions in partially ionized plasma 

Keeping the quasi-neutrality limit, we now discuss the I A wave damping in plasmas 
comprising neutrals as well. In view of the results presented above, the electron-ion 
friction terms in Eqs. ([T]) and d2J) will cancel each other out and in a few steps one derives 
the following dispersion equation containing the collisions of plasma species with neutrals 
and vice versa: 




(6) 
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Figure 1: Frequency u r and absolute value of the IA mode damping |7| in terms of the 
number density of neutrals. Details of the mode behavior in the region A are better seen 
in the linear scale (small figure inside). 

In the derivation, the ion and neutral thermal terms are neglected. The ion thermal terms 
would give the modified mode frequency u> 2 = k 2 c 2 s (l + Ti/T e ). Hence, even if T e = T{ the 
wave frequency is only modified by a factor 2 1//2 . The neutral thermal terms are discussed 
further in the text. Note that in deriving Eq. ((H]), the momentum conservation condition 
v ie = m e v ei /mi is nowhere used: the e-i and i-e friction terms exactly vanish in view of 
Eq. ©. 

Equation ([8]) is solved numerically for a plasma containing electrons, protons, and 
neutral hydrogen atoms using the following set of parameters: T e = 4 eV, n$ = 10 18 m -3 , 
k — 10 m _1 , with [10] a en = 1.14- 10 -19 m -2 . The neutral density is varying in the interval 
10 16 — 10 23 m -3 . The ion and hydrogen temperatures are taken T» =T n = T e /20, satisfying 
the condition of their small thermal effects. This also gives [TTJ, Oi n = 2.24 ■ 10~ 18 m -2 . 
The results are presented in Fig. 1. The IA mode propagates in two distinct regions A 
and B. 

Only a limited left part of the region A would correspond to the 'standard' IA wave 
behavior in a collisional plasma: the mode is damped and the damping is proportional 
to the neutral number density. Hence, in this region it may be more or less appropriate 
to use the approximate expressions for the friction force, like (in the case of electrons) 
F e ~ m e nov en v e i. However, this domain is very limited because in the rest of the domain 
the frequency drops and the mode becomes non-propagating for n n Q > 3.8 ■ 10 19 m -3 (this 
is the lower limit of the region C in Fig. 1). 
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Figure 2: Frequency u r and absolute value of the IA mode damping |7| in terms of the 
wave number. The line kc s shows part of the graph of the ideal mode. Details of the 
domain a are better seen in the linear scale (small figure inside). 

Increasing the neutral number density, after some critical value (in the present case this 
is around n„o — 10 20 m~ 3 ) the IA mode reappears again in the region B, with a frequency 
starting from zero. For even larger neutrals number densities, the mode damping in fact 
vanishes completely and the wave propagates freely but with a frequency that is many 
orders of magnitude below the ideal case kc s ~ 196 kHz. This behavior can be explained 
in the following manner. For a relatively small number of collisions the IA mode is weakly 
damped because initially neutrals do not participate in the wave motion and do not share 
the same momentum. Increasing the number of neutrals, the damping may become so 
strong that the wave becomes evanescent. However, for much larger collision frequencies 
(i.e., for a lower ionization ratio), the tiny population of electrons and ions is still capable 
of dragging neutrals along and all three components move together. The plasma and the 
neutrals become so strongly coupled that the two essentially different fluids participate 
in the electrostatic wave together. In this regime, the stronger the collisions are, the less 
wave damping there is! Yet, this a bit counter-intuitive behavior comes with a price: the 
wave frequency and the wave energy flux becomes reduced by several orders of magnitude. 

Similar effects may be expected by varying the wave-length. The previous role of the 
varying density of neutrals is now replaced by the the ratio of the mean free path of a 
species Xfj = v T j/uj (with respect to their collision with neutrals) and the wavelength. 
This ratio now determines the coupling between the plasma and the neutrals. The mode 
behavior is directly numerically checked by fixing n n o = 10 20 m -3 , no = 10 18 m -3 , and 
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Figure 3: The two lines give the lower and upper values of the neutrals' density n n0 
between which, for the given plasma density n , the IA mode does not propagate. 

for other parameters same as above. For these parameters we have A/ e = v Te /u en = 0.09 
m, and Xfi — v T i/u in = 0.004. The numerical results are presented in Fig. 2 for k varying 
in the interval 0.2 — 80 m _1 . The mode vanishes in the interval c, between k ~ 10 m _1 
and k ~ 25.6 m -1 . The explanation is similar as before. Note that for k = 0.2 m -1 (in 
the region a) we have ui r ~ 390 Hz, and this is about one order below kc s . Compared to 
the mode behavior in Fig. 1, this implies that the mode in the present domain a is in the 
regime equivalent to the domain B from Fig. 1; here, in Fig. 2, these large wave-lengths 
imply well coupled plasma-neutrals, where the frequency is reduced and the damping is 
small. The region a is also given separately in linear scale together with the dotted line 
describing the ideal mode kc s . Clearly, in general the realistic behavior of the wave is 
beyond recognition and completely different as compared to the ideal case. 

After checking for various sets of plasma densities, it appears that the evanescence 
region reduces and vanishes for larger plasma densities no- This is presented in Fig. 3 for 
the same parameters as above, by taking k = 10 m _1 , but for a varying plasma density 
n . The two lines represent boundary values of the number densities of neutrals, for the 
given plasma density, at which the IA mode vanishes; for the neutrals densities between 
the two lines the IA mode does not propagate. The symbols * on the two lines denote 
the boundaries of the region C from Fig. 1. It is seen that for the given case the IA 
mode propagates without evanescence for the plasma densities above n = 3.8 ■ 10 18 m~ 3 . 
Physical reason for a larger non-propagating domain for low plasma density is obvious, 
namely the tiny plasma population is less efficient in inducing a synchronous motion of 
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Figure 4: Values of the wave-number, in terms of the plasma density, for which the IA 
wave becomes evanescent. In the region between the lines the mode does not propagate. 

neutrals. In the other limit, the opposite happens and the forbidden region eventually 
vanishes. 

A similar check is done by varying the wave-number and the plasma density, and the 
result is presented in Fig. 4 for a fixed n„o = 10 20 m -3 . The lines represent the values 
(no, k) at which the IA wave becomes evanescent. There can be no wave in the region 
between the lines. On the other hand, there is no evanescence for the plasma density 
above Uq = 1.2 ■ 10 19 m -3 . Here * denote the boundaries of the region c from Fig. 2. 

All these results clearly indicate that in practical measurements in laboratory and 
space plasmas, the IA mode can hardly be detected and recognized as the IA mode unless 
collisions are correctly taken into account (using full friction terms), and the mode is 
sought in the corresponding domain which follows from our Eq. flH]). 

2.3 Thermal effects of neutrals 

Keeping the pressure terms for ions and neutrals yields the following dispersion equation 



Here, v 2 = c 2 + v \ { . Without collisions, this yields two independent modes, viz. the 
ion-acoustic mode and the gas thermal (GT) mode, [yj 2 — k 2 v 2 n ) (cu 2 — k 2 v 2 s ) = 0. The 
collisions couple the two modes, and in order to compare with the previous cases we solve 
Eq. (ED for k = 10 m -1 , n = 10 18 m~ 3 , T e = 4 eV, % = T e /20, and in terms of the 





(9) 
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Figure 5: The real part of the frequency of damped gas thermal mode in terms of the 
number density of neutrals and for several temperatures of the neutrals gas. 

density and temperature of neutrals. For a low thermal contribution of neutrals (i.e., a 
low neutral temperature, or/and heavy neutral atoms) the previous results remain valid. 
Larger values of v Tn introduce new effects, this is checked by varying the temperature T n . 
The ion thermal terms do not make much difference, as explained earlier. The real part of 
the frequency u g of the gas thermal mode is presented in Fig. |5j and this only in a limited 
region that includes the evanescence area C from Fig. [TJ The damping is not presented 
but the mode is in fact heavily damped. 

The explanation of the figure is as follows. The starting solution for T n = is in fact 
the line u g = 0, and this case would correspond to the the IA mode from Fig. [TJ For 
some finite T n there appears the GT mode. For a low gas temperature the mode becomes 
evanescent for a higher density of neutrals (the dot and dash-dot lines in Fig. [5]). This 
evanescence is accompanied with the previously discussed evanescence and re-appearance 
of the IA mode (described earlier and no need to be presented here again). However, 
for still larger T n , the IA and GT modes become indistinguishable and propagate as one 
single mode. This is presented by the two upper (the full and dashed) lines in Fig. [5j 
that go up for large enough n n Q. Also given are the corresponding ideal values kv Tn that 
appear to be much above the actual wave frequency u g in such a collisional plasma, but 
this remains so only until the neutral density n n o exceeds some critical value. After that 
the wave in fact behaves as less and less collisional and the wave frequency is increased. 
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3 IA wave instability in inhomogeneous partially ion- 
ized plasma 

3.1 Fluid description in collisional plasma 

In the previous text, collisions were shown to yield damping of the IA mode. However, if 
the plasma is inhomogeneous, implying the presence of source of free energy in the system, 
a drift-type instability of the IA wave may develop if there is a magnetic field B = B e z 
present, and the electrons (ions) are magnetized (un- magnetized). The magnetic field 
introduces a difference in the parallel and perpendicular dynamics of the magnetized 
species so that the continuity condition in this case can be written as 



driji 



+ 71,-0 V " Vji + Vji ■ Vrij-o = 0. (10) 



Here, V = Vj_ + V 2 . For the un-magnetized species the direction of the wave plays 
no role so that V — > ik, k 2 = k 2 + k 2 . On the other hand, for the equilibrium gradi- 
ent along the x-axis and for perturbations of the form ~ f{x) exp(— iut + ik y y + ik z z), 
where \(df /dx)/ f\, \ (drijo/dx)/no\ -C k y , we apply a local approximation, and for ions the 
last term in Eq. ffTUl) vanishes because of the assumed geometry. The ions' dynamics is 
basically the same as in the previous sections. 

The electron momentum equation (j2J) will now include the Lorentz force term —en e v e x 
B. Repeating the derivation from Ref. [3], the total perpendicular electron velocity can 
be written as 



v e ± 



1 + i&aVfi* 



1 _ v en aV ± (j) v 2 T u en aV^n e w| V ± n e 

— e z X V_L0 + — TT e z x 



Til 



B fi e Bq VL 2 e n e £l e n e 

In the direction along the magnetic field vector, the perturbed electron velocity is 

v ezl = ih< U2 + UL ( ^ - (12) 



Here, a — uj (uj + iu ne ), and for magnetized electrons, |i/g n a 2 /f2g| <C 1 in the denominator 
in Eq. ( fiTT) . Using these equations in the continuity condition (TTDjl for electrons one 
obtains 

n el + iD p + iD z [bj 2 + h>l e )/{uj 2 - iv ne bS) efa 



n lo + iD p + iD z {u 2 + v 2 e )/(u 2 - iv ne uj) nT e 
D p = v en ak 2 p 2 e , D z = k 2 vlJu en , p e = v Te /Vt e . 



(13) 
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The term D p describes the effects of collisions on the electron perpendicular dynamics 
and is usually omitted in the literature. However, as shown in a recent study j3], it can 
strongly modify the growth rate of the drift and IA-drift wave instability in the limit of 
small parallel wave- number k z . 

Neglecting the neutral dynamics is equivalent to setting v ne = 0. This yields a = 1, 
and Eq. (fT5|) becomes identical to the corresponding expression in Refs. [21 12] • For a 
negligible D p , Eq. (1T5|) becomes the same as the corresponding equation from Ref. [T5] . 
For negligible ion thermal effects, the final dispersion equation reads 

k 2 c 2 s = + iD p + iD z (u 2 + v 2 e )/(u 2 - iv ne uj) ,^ 
tu 2 tu + iD p + iD z (u 2 + i> 2 e )/(uj 2 - iu ne u) 

Equation (I14p can be solved numerically keeping in mind a number of conditions used 
in their derivations, like smallness of the plasma beta to remain in electrostatic limit, 
smallness of the parallel phase velocity as compared to the electron thermal speed because 
of the massless electrons limit, also the ratio D p /D z should be kept not too big or too 
small in order to have the assumed effects of electron collisions in perpendicular direction. 
We plan to compare this collisional instability with the kinetic instability due to the 
presence of the density gradient. Therefore, the wave frequency should be below the 
electron diamagnetic frequency etc. 

We solve Eq. (j!4p for an electron-argon plasma in the presence of parental argon atoms. 
As an example we take T e = 4 eV, T { = T n = T e /30, n = 10 15 m" 3 , B = 1.2 ■ 1(T 2 T, 
k = 500 m -1 , L n = 0.05 m, and take several values for the density of neutrals. The result 
in terms of the angle of the propagation 9 = aicta.n(k z /k y ) is presented in Fig. 5. The 
three lines (full for the real part of the frequency, and dashed for the growth rates) are 
for n„o = 10 19 , 10 18 , 10 17 m -3 . It is seen that i) the instability is angle dependent and 
there exists an angle of preference and an instability window in terms of 9 within which 
the mode is most easily excited, ii) this angle of preference is shifted towards smaller 
values for lower values of the neutral density, and iii) in the same time the instability 
window becomes considerably reduced. This shows an interesting possibility of launching 
the IA-drift wave in a certain direction by simply varying the pressure of the neutral gas. 

Varying the density scale length L n = (dno/dx)' 1 the wave frequency may become 
above u* e and in this case the instability vanishes. As an example, this is demonstrated 
in Fig. [7] for the parameters corresponding to the line II from Fig. |6] and for the angle 9 
at the maximum growth rate. The growth rate changes the sign for u ~ u* e . 
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Figure 6: Real part of the frequency from Eq. (1 1 4jl (full lines) and the corresponding 
growth rates (dashed lines), both normalized to the electron diamagnetic drift frequency, 
for three values of neutral number density. The lines I, II, III correspond (respectively) 
to n n Q = 10 19 , 10 18 , 10 17 m~ 3 . The line 7^ is the kinetic growth-rate from Eq. (TT9|) (for 
the same parameters as line II). 
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Figure 7: The real and imaginary parts of the frequency for the line II from Fig. [6j in 
terms of the characteristic density inhomogeneity scale length L n = (dn /dx)~ l , and for 
the angle 9 at the maximum on Fig. EJ 
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3.2 Comparison with collision-less kinetic gradient-driven IA 
wave instability 

Keeping the same model of magnetized (un- magnetized) electrons (ions), within the ki- 
netic theory the perturbed number density for electrons can be written as 

2 1 



n 



e<pi 



7T 
I I - 



U) — CJ* 



exp 



-co 2 /(2k 2 z v 



2. .2 > 

Te/ 



(15) 



In the derivation of Eq. ffl5l) the electron Larmor radius corrections are neglected in terms 
of the type I n {b) exp(— b), b = k\p 2 , where I n denotes the modified Bessel function of the 
first kind, order n, and only n = terms are kept for the present case of frequencies much 
below the gyro-frequency 

The ion number density can be calculated using the kinetic description for un-magnetized 
species, the derivation is straight-forward and it yields [15] 



Tig 

n i0 



e<pi 



rriivli 



(16) 



Here, J(rj) = [r]/(2n) 1 ^ 2 ] J c d(exp(—( 2 /2)/(r] — £) is the plasma dispersion function, and 
( = v/v Ti . In the case \r]\ 3> 1, and assuming |i?e(?7)| 3> Im(r]), an expansion is used for 
J {if). This together with the quasi-neutrality yields the kinetic dispersion equation for 
the IA-drift wave: 



L2 2 

A(u, k) = 1 - 



i(n/2) 



1/2 



td - CO* 



k z Vn 



exp 



-tu 2 /(2k 2 v 



2„.2 > 
Te/ 



+ 



J. I" Jv 1)t". 



The real part of Eq. (1171) yields the spectrum 



k 2 c 2 



1 + (1 + Y2T.IT,) 



exp 



1/2' 



-u 2 /{2k 2 v. 



2_,2 > 
TiJ 



(17) 



;is) 



The kinetic growth rate is given by 

7^ ~ —ImA/{dReA/duj 



(n/2Y/ 2 ul 



2k 2 c 2 



x 



k z Vn 



exp 



2_.2 \ 

Te/ 



+ ThT exp 



-u\l{2k 2 v, 



2„.2 > 

TiJ 



(19) 



Here, the index k is used to denote kinetic expressions. The electron contribution in 
Eq. ffl9|) yields a kinetic instability provided that ujk < w* e - 

Equation ( TT9|) is solved numerically and compared with the growth rate obtained from 
the collisional IA-drift mode (jHJ). For a fixed k = 500 m _1 as in Figs. [6] and [TJ the 
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normalized frequency Uj-foj^ = 0.485, and the result for the growth rate is presented 
by the line 7^ in Fig. [6] for the parameters corresponding to the line II from the fluid 
analysis (i.e., for n n Q = 10 18 m~ 3 ). The larger kinetic growth rate appears also to be 
angle dependent, yet with a much wider instability window as compared to the collisional 
gradient driven instability obtained from the fluid theory. 

4 Summary 

The analysis of the ion acoustic wave presented here shows the importance of collisions in 
describing the wave behavior. Without a proper analytical description, the identification 
of the mode in the laboratory and space observations may be rather difficult because one 
might fruitlessly search for the wave in a very inappropriate domain, as can be concluded 
from the graphs presented here, and in particular from Fig. 2. Not only the wave frequency 
may become orders of magnitude below an expected ideal value, but also the mode may 
completely vanish. A similar analysis of the effects of collisions may be performed for 
other plasma modes as well, like the Alfven wave etc, as predicted long ago in classic 
Ref. [16J. The impression is that these effects are frequently overlooked in the literature, 
hence the necessity for the quantitative analysis given in the present work that can be 
used as a good starting point for an eventual experimental check of the wave behavior 
in collisional plasmas. Particularly interesting for experimental investigations may be 
the angle dependent mode behavior given in Sec. 3, where it is shown that the strongly 
growing mode may be expected within a given narrow instability window in terms of the 
angle of propagation. Comparison with the kinetic theory shows a less pronounced angle 
dependent peak, yet this kinetic effect can effectively be smeared out in the presence 
of numerous collisions, that are known to reduce kinetic effects in any case, and the 
sharp angle dependence that follow from pure fluid effects should become experimentally 
detectable. 
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